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ABSTRACT 


The  basic  equations  of  motion  for  a  point  mass  under  the  influence 
of  gravity,  aerodynamic  drag,  end  rocket  forces  are  used  to  simulate  a 
three-dimensional,  three-degree  of  freedom,  helicopter  emergency  egress 
trajectory.  The  resulting  equations  of  motion  were  solved  by  a  fourth 
order  Runge  Kutta  method,  programmed  in  Fortran  for  the  Univac  Solid 
State  90  digital  computer. 

Included  in  this  report  is  the  computer  program  which  was  used 
to  obtain  a  numerical  solution  to  the  equations  of  motion,  the  solu¬ 
tion  to  three  test  cases,  and  a  comparison  with  an  experimental  tra^ 
jectory . 
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INTRODUCTION 


Operational  experience  with  Army  helicopters  reveals  many  in¬ 
stances  where  Ln-flight  hazards  to  crew  personnel  may  be  effectively 
resolved,  with  a  resultant  increase  in  flight  safety,  through  the  use 
of  an  emergency  escape  system. 

In  this  report,  the  basic  equations  of  motion  for  a  point  mass 
under  the  influence  of  gravity,  aerodynamic  drag,  and  rocket  forces 
are  used  to  simulate  a  three-dimensional,  three-degree  of  freedom, 
helicopter  emergency  egress  trajectory.  The  resulting  equations  of 
motion  were  solved  by  a  fourth  order  Runge  Kutta  method,  programmed 
in  Fortran  for  the  Univac  Solid  State  90  digital  computer.  The  com¬ 
puter  program  which  was  used  to  obtain  a  numerical  solution  to  the 
equations  of  motion  is  presented  in  Appendix  A. 

The  purpose  of  this  report  is  to  document  the  analytical  capa¬ 
bility  necessary  to  evaluate  the  egress  trajectories  associated  with 
such  emergency  escape  systems. 


ASSUMPTIONS 


1.  The  seat-man  combination  is  treated  as  a  point  mass. 

2.  Time,  t  =  0,  begins  at  the  end  of  the  catapult  stroke,  at 
which  time  the  catapult  has  imparted  some  initial  velocity  to  the 
seat-man  combination. 

3.  Aerodynamic  drag  is  proportional  to  the  n^  power  of  the 
velocity.  In  the  test  cases  presented,  n  is  2. 

4.  Rocket  thrust  is  constant  in  both  magnitude  and  direction 
throughout  the  burn  time. 

5.  The  helicopter  pitch  angle  is  zero  at  the  time  of  ejection. 

6.  Parachute  forces  are  omitted. 
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DISCUSSION 


•/ 


The  helicopter  emergency  egress  trajectory  analyzed  in  this  re¬ 
port  is  accomplished  initially  by  a  catapult  (which  ejects  the  seat-man 
laterally)  and,  after  a  suitable  delay,  by  a  rocket  (which  raises  the 
escapee  to  a  height  sufficient  for  safe  recovery  by  parachute)  .  The 
analysis  begins  at  the  end  of  catapult  stroke  (t  =  0) ,  at  which  time 
the  seat-man  has  acquired  an  initial  velocity  due  to  both  the  cata¬ 
pult  and  the  velocity  of  the  supporting  vehicle.  The  magnitude  and 
direction  of  these  velocities  are  inputs  to  the  program. 

The  trajectory  can  be  thought  of  in  terms  of  the  following  three 
stages : 


1.  Initial  "Free"  Flight  Stage  -  The  motion  of  the  seat -man, 
determined  by  the  initial  conditions,  continues  under  the  influences 
of  only  gravity  and  drag  up  to  time  t^. 

2.  Rocket  Stage  -  At  time  tj_,  a  rocket  force  (R)  becomes  effec¬ 
tive.  This  force,  constant  in  time  and  direction,  terminates  at  time 
ti  +  tb»  The  delay  between  ejection  (t  =  0)  and  rocket  initiation  (tj) 
allows  the  seat-man  sufficient  time  to  clear  the  rotor  blades  before 
being  propelled  upward.  The  magnitude  and  direction  of  R,  t^,  and  t^ 
are  all  inputs  to  the  program. 

3.  Final  "Free"  Flight  Stage  -  At  time  t^  +  tb,  the  seat-man 
continues  its  trajectory,  influenced  only  by  gravity  and  drag. 

The  trajectory  is  terminated  automatically  in  either  of  two  ways: 
(1)  when  ground  impact  is  achieved,  i.e.,  when  the  vertical  coordinate, 
X3,  first  becomes  negative;  or  (2)  when  a  terminal  time  read  into  the 
computer  as  input  data  is  reached;  whichever  occurs  first. 


THE  COORDINATE  SYSTEM 


A  right-handed  orthogonal  coordinate  system  (Figure  1)  is  employed 
with  the  mutually  orthognal  axes  being  defined  in  the  following  manner. 

The  +X3  axis  is  chosen  to  be  in  a  direction  antiparallel  to  that 
of  the  acceleration  due  to  gravity.  The  roll  axis  of  the  helicopter 
was  chosen  as  the  x^  axis,  with  the  fore  end  of  the  helicopter  de¬ 
fining  the  +  direction.  The  x 2  axis  is  then  so  oriented  as  to  com¬ 
plete  the  right-handed  orthogonal  system.  The  xj_,  X£,  and  X3  cordinate 
axes  are  referred  to  as  x,  y,  and  z,  respectively,  in  the  computer 
input/output  formats. 
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Figure  1.  The  Coordinate  System 


All  directional  inputs  to  the  program  are  read  into  the  computer 
in  terms  of  the  angles  6  and  <p  (the  colatitude  and  azimuth,  respec¬ 
tively)  referenced  with  respect  to  this  space  set  of  axes  fixed  at  the 
instant  of  ejection. 


i 
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EQUATIONS  OF  MOTION 

Summing  forces  (see  Appendix  B  for  symbols  used) 

£  Fi^  =  mxi 


j 


where 


Fi 

Fi 


and 


*Dots  represent  differentiation  with  respect  to  time. 


(l)* 


(1)  =  X-mg  for  i  -  3 

l  0  for  i  -  1,  2 

(Gravity) 

(2)  =  -k  vn  vj/v  =  -k  vn_l  V£ 

(Drag) 

(3)  =  for  tL  *  t  *  tx  +  tb 

(.0  otherwise 

(Rocket) 

«<  • 
xi  =  vi 

V  =  (vi2  + 

V22  +  V32) 2 

(2) 

(3) 
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Combining  Equations  1,  2,  and  3  for  n  =  2, 


mv£ 


k  v  +  F^3^ 


(4) 


METHOD  OF  SOLUTION 


Equation  4  represents  the  three  simultaneous  differential  equa 
tions : 

*1  “  "  £  *  V1  +  Fl(3> 

$2  =  '  m  v  v2  +  F2(3) 

V3  =  '  g  -  i  v  V3  +  F3(3) 

These  three  equations,  when  coupled  with 


*1  =  V1 
x2  =  v2 

x3  =  v3 

result  in  six  first  order  differential  equations  which  are  integrated 
numerically  using  a  fourth  order  Runge  Kutta  method  to  obtain  the 
velocities  and  coordinates  as  functions  of  time: 

v±  =  v^t) 
xi  =  xi(fc) • 


COMPUTER  INPUT  FORMAT 


The  data  input  format  for  the  computer  program  is  listed  in 
Table  I.  The  input  data  are  read  into  the  computer  on  a  series  of 
18  cards,  in  the  order  shown.  The  first  card  titles  each  particular 
trajectory  simulation.  This  information  is  punched  into  the  first 
portion  of  the  card,  beginning  in  column  1.  If  a  title  is  not  de¬ 
sired,  a  blank  card  must  be  substituted.  The  next  17  cards  contain 


4 


the  necessary  input  data.  Parameter  titles  are  punched  into  columns 
1  through  45.  Columns  47  through  90  are  reserved  for  the  numerical 
values  of  the  input,  column  46  containing  the  input  sign. 

An  example  of  a  suggested  format  is  shown  in  the  initial  portion 
of  the  computer  printout  (Figure  2)  for  Test  No.  1.  The  Fortran  symbols 
designating  the  various  inputs  to  the  computer  program  are  also  listed 
in  Table  I. 


TABLE  I.  Input  Data 

Test  No.  _ 


Fortran  Symbols 

Parameter 

Unit 

D(l) 

Initial  height 

ft 

D(2) 

Vehicle  velocity 

ft/sec 

D(3) 

Catapult  velocity 

ft/sec 

D<4) 

Rocket  thrust 

lb 

D(5) 

Vehicle  colatitude,  0 

degree 

D(6) 

Vehicle  azimuth,  <p 

degree 

D(7) 

Catapult  colatitude,  $ 

degree 

D(8) 

Catapult  azimuth,  <p 

degree 

D(9) 

Rocket  colatitude,  0 

degree 

D(10) 

Rocket  azimuth, 

degree 

D(U) 

Rocket  initiation  time 

sec 

D(12) 

Rocket  burn  time 

sec 

D(13) 

Ejected  mass 

slugs 

D(14) 

Drag  coefficient 

lb(sec/ft)n 

D(1 5) 

Drag  exponent 

dimensionless 

D(16) 

Time  increment 

sec 

D(17) 

Terminal  time 

sec 

Explanation  of  the  input  data  follows . 

D(l)  -  The  initial  h  ight  of  the  seat-man  configuration. 

D(2)  -  The  magnitude  of  the  helicopter  velocity  at  t  =  0. 

D(3)  -  Magnitude  of  the  velocity  imparted  to  the  seat-man  by  the  cata¬ 
pult  at  t  =0. 

D(4)  -  Rocket  thrust. 

D(5)T  The  direction  of  the  vehicle  velocity  at  t  =  0  referenced  with 
D(6)J  respect  to  the  space  set  of  axes. 

D(7)]  The  direction  of  the  catapult  velocity  referenced  with  respect. 
D(8)J  to  the  space  set  of  axes. 
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HELICOPTER  EGRESS  TRAJECTORY  ANALYSIS 
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Figure  2.  Initial  Portion  of  Computer  Printout  for  Test 


D(9)  1  Direction  of  rocket  thrust  referenced  with  respect  to  the 
D(10)J  space  set  of  axes. 

D(ll)  -  Rocket  initiation  time. 


D(12)  -  Rocket  burn  time. 

D (13)  -  Mass  of  seat-man  combination, 

D(14)  -  Drag  coefficient  equal  to  pAC])/2,  where  p  is  the  air  density 
(slugs/ft-^);  A,  the  seat  area  (ft^)  ;  and  Cd,  an  aerodynamic 
coefficient  (dimensionless)  taken  equal  to  1. 

D(15)  -  Velocity  exponent  in  the  velocity  dependent  drag  force. 

D(16)  -  Compute  and  print  interval. 

D(17)  -  Time  at  which  the  program  terminates  if  ground  impact  has  not 
occurred . 


RESULTS  OF  COMPUTER  ANALYSIS 


Approximately  twelve  computer  tests  were  run;  the  results  of  three 
(tests  1,  5,  and  6)  are  given  in  Figures  3,  4,  and  5.  The  input  for 
these  tests  is  presented  in  Table  II. 

The  computer  simulation  was  compared  to  data  from  experimental 
trajectories.*  Test  1  is  such  a  comparison.  The  input  values  were 
chosen  to  correspond  to  those  of  the  actual  ejection.  Figure  2  is 
the  initial  portion  of  the  computer  printout  for  this  test.  Coor¬ 
dinates  (x,  y,  z)  are  expressed  in  feet;  velocities  (vx,  Vy,  vz)  ,  in 
ft/sec;  and  accelerations  (Ax,  Ay,  Az) ,  in  ft/sec^.  The  trajectory 
is  referenced  with  respect  to  the  point  of  ejection.  Comparison  be¬ 
tween  experimental  and  computer  results  is  presented  in  Table  III. 

Tests  5  and  6  were  run  using  the  same  aerodynamic  and  rocket 
data,  but  with  different  initial  heights  and  vehicle  velocities. 

A  suitable  method  of  graphically  portraying  the  three-dimensional 
trajectory  is  shown  in  Figures  3,  4,  and  5.  Two  graphs  (A  and  B)  de¬ 
termine  each  trajectory.  The  first  (A)  is  a  graph  of  altitude  (z)  vs 
time  (t) ,  with  the  various  events  indicated  along  the  time  axis.  The 
second  (B)  is  a  graph  of  the  horizontal  components  x  vs  y,  again  with 
the  various  events  indicated.  The  two  graphs  are  sufficient  to  de¬ 
termine  the  complete  space-time  trajectory. 

*  M.  Weinstock,  "Emergency  Escape  Systems  for  Army  Helicopters," 
Aerospace  Medicine,  Vol  36,  No.  3,  March  1963. 
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TABLE  II.  Input  for  Tests  I,  5,  and  6 


Parameter 

Test  1 

Test  5 

Test  6 

Initial  height 

(ft) 

5.0 

50.0 

50.0 

Vehicle  velocity 

(ft/sec) 

0.0 

50.0 

50.0 

Catapult  velocity 

(ft/sec) 

60.0 

60.0 

60.0 

Rocket  thrust 

(lb) 

3400 .0 

3400.0 

3400.0 

Vehicle  colatitude 

(degree) 

90.0 

180.0 

0.0 

Vehicle  azimuth 

(degree) 

0.0 

0.0 

0.0 

Catapult  colatitude  (degree) 

90.0 

90.0 

90.0 

Catapult  azimuth 

(degree) 

90.0 

90.0 

90.0 

Rocket  colatitude 

(degree) 

36.25 

36.25 

36.25 

Rocket  azimuth 

(degree) 

0.0 

0.0 

0.0 

Rocket  initiation 

time  (sec) 

0.3 

0.3 

0.3 

Rocket  burn  time 

(sec) 

0.35 

0.35 

0.35 

Ejected  mass 

(slugs) 

9.3243 

9.3243 

9.3243 

Drag  coefficient 

(lb/sec/ft)n 

0.006 

0.006 

0.006 

Drag  exponent  (dimensionless) 

2.0 

2.0 

2.0 

Time  increment 

(sec) 

0.05 

0.05 

0.05 

Terminal  time 

(sec) 

10.0 

10.0 

10.0 

TABLE 

III.  Comparison  between  Experimental 

and  Computed  Trajectory 

Max  Height 

Max  vx 

Max  Vy 

Max  vz 

(ft) 

(ft/sec) 

(ft/sec) 

(ft/sec) 

Experimental 

120  to  130 

72 

60 

100 

Computed 

119 

78 

60 

84 

No  horizontal  displacement  data  are  available  for  comparison, 
and  no  meaningful  comparison  between  acceleration  levels  is  possible 
because  of  the  following  assumptions  made  in  the  analysis: 

1.  The  catapult  which,  in  the  experimental  trajectory,  caused 
the  maximum  acceleration  is  assumed  to  impart  an  initial  velocity  only; 

2.  Rocket  thrust  is  assumed  constant  over  the  burn  time. 
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CONCLUSIONS 


Comparison  with  experimental  data  indicates  that  the  analysis 
presented  in  this  report  adequately  simulates  experimental  helicopter 
egress  trajectories. 


RECOMMENDATIONS 


It  is  recommended  that 

1.  The  analytical  capability  provided  here  be  utilized  in 
future  analyses  of  helicopter  egress  trajectories  and  any  other  tra¬ 
jectory  simulation  fulfilling  the  stated  assumptions. 

2.  This  analysis  be  extended  to  include  the  effect  of  para¬ 
chute  forces  and  rotation  in  the  pitch  plane. 
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APPENDIX  A 
THE  COMPUTER  PROGRAM 


C  HELICOPTER  EGRESS  TRAJECTORY  ANALYSIS 

DIMENSION  A(10),B(10),C(10),D(25),F(10),HED(9),S(10),U(10) 
PRINT  402, 

PRINT  401, 

READ  202,HED 
PRINT  203, HED 
PRINT  401, 

DO  300  I#l,17 
READ  200, HED, D(  I) 

300  PRINT  204, HED, D( I) 

PRINT  401, 

PRINT  400, 

C  CONSTANTS 

T2#D(11)+D(12) 

D(9)#D(9)-90. 

D(10)#D(10)-90. 

DO  60  I#5,10 
60  D(I)#.017453*D(I) 

DO  70  I#5 ,8 
70  A( I )#SIN(D( I ) ) 

DO  80  I#5,8 
80  B(I)#COS(D(I)) 

C1#SIN(D(7 )+D(9 ) ) 

C2#C0S(D(7)+D(9) ) 

C3#SIN(D(8)+D(10)) 

C4#C0S(D( 8 )+D(10 ) ) 

C5//D(2)*A(5)-'B(6) 

C6#D(2)*A(5)*A(6) 

C7#D(2 )*B( 5 ) 

C8#D(3)*A(7)*B(8) 

C9#D ( 3 )*A ( 7 )*A ( 8 ) 

C10#D(3)*B(7) 

C11#D(4  )*C1*C4/D(13 ) 

C12#D(4)*C1*C3/D(13) 

C13#D(4)*C2/D(13) 

C14#-D(14)/D(13) 

C15#2.*D(16) 

C16#.5*D(16) 

C  INITIAL  CONDITIONS 
J#1 
T#0. 

XD0C5+C8 

YD#C6tC9 

ZD0C7+C1O 

X#0. 

Y#0  • 

Z#D(1) 

F(1)#XD 
F( 2  )#YD 
F(3)#ZD 
F(4)#X 
F(5)#Y 

r(6)#z 


16 


C  ENTER  EVALUATION  ROUTINE 
100  XDJSfF(l) 

YD#F(2) 

ZD#F( 3) 

X#F(4) 

Y#F(5) 

Z#F(6) 

V#SQRT(XD*XD+YD*YD+ZD*ZD) 

IF(T-D(11))2,1,1 

1  IF(T-T2)3,3,2 

2  RX#0. 

RY#0. 

RZ#0. 

GO  TO  4 

3  RX0C11 
RY#C12 
RZ0C13 

4  VC14#C14*V**(D(15)-1.) 

S(  1  )# VC14*XD+RX 
S(2)#VC14*YD+RY 

S(  3  )#VC14*ZD+RZ-32 .174 
S(4)#XD 
S(5)#YD 
’ S(6 )#ZD 
F(1)#S(4) 

F(2)#S(5) 

F(3)#S(6) 

F(4)#X 

F(5)#Y 

F(6)#Z 

IF(T)6,5,6 

5  PRINT  201,T,F(4),F(5),F(6),F(1),F(2),F(3),S(1),S(2),S(3) 

C  EXIT  EVALUATION  ROUTINE 

6  GO  TO  (10, 20,30,40, 50), J 

C  ENTER  RUNGE  KUTTA  ROUTINE 

10  DO  11  I#l,6 

ummn 

C( I )#D(16 )*S( I ) 

11  F( I)#U(  I )+. 5*C(I ) 

T0T+C1S 

J#2 

GO  TO  100 

20  DO  21  1>U,6 
C(I)#C(I)+C15’-S(I) 

21  F(I)#U(I)+C16*ba) 

m 

GO  TO  100 

30  DO  31  I#l,6 

C( I )#C( I )+C15*S( I ) 

31  F(I)#U(I)+D(16)*S(I) 

T0T+C16 

m 

GO  TO  100 
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40  DO  41  I#lf6 

41  F(I)#U(I)+(C(I)+D(16)*S(I))/6. 

J#5 

GO  TO  100 

EXIT  RUNGE  KUTTA  ROUTINE 

50  PRINT  201,T,F(4),F(5),F(6),F(1),F(2),F(3),S(1),S(2),S(3) 

Ml 

IF(Z)102,102 ,101 

101  IF(T-D(17) >100,102,102 

102  STOP 

200  FORMAT (9A5,  F16.5) 

201  F0RMAT(10F13.2) 

202  FORMAT (9A5 ) 

203  FORMAT( 59X.9A5) 

204  FORMAT (34X, 9 A5,  F16.5) 

400  FORMAT ( 10X , 1HT , 12X , 1HX , 12X , 1HY , 1?X , 1HZ , 12X , 2HVX , 11X , 2HV  Y , 11X , 2HVZ , 
111X,2HAX,11X,2HAY,10X,2HAZ) 

401  FORMAT( 2/ ) 

402  F0RMAT(46Xl 37HHELIC0PTER  EGRESS  TRAJECTORY  ANALYSIS) 

END 


APPENDIX  B 


SYMBOLS  USED  IN  DERIVATION  OF  EQUATIONS  OF  MOTION 


Symbol 

Definition 

Unit 

F 

Force 

lb 

g 

Acceleration  due  to  gravity 

f t/sec^ 

i 

Subscript  having  values  1,  2,  3, 
corresponding  to  the  three  mutually 
perpendicular  coordinate  directions 

dimensionless 

j 

Superscript  having  values  1,  2,  3, 
referring  to  the  various  force  con¬ 
tributions 

dimensionless 

k 

Drag  coefficient 

lb(sec/ft)n 

m 

Mass  of  seat-man 

slugs 

n 

Drag  exponent 

dimensionless 

R 

Rocket  thrust 

lb 

t 

Time 

sec 

tl 

Rocket  initiation  time 

sec 

tb 

Rocket  burn  time 

sec 

V 

Velocity 

ft/ sec 

X 

Displacement 

ft 

iy 
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